Detection of Massive Tidal Tails around the Globular Cluster 
Pal 5 with SDSS Commissioning Data 
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ABSTRACT 

We report the discovery of two well-defined tidal tails emerging from the sparse remote glob- 
ular cluster Palomar 5. These tails stretch out symmetrically to both sides of the cluster in the 
direction of constant Galactic latitude and subtend an angle of 2.6° on the sky. The tails have 
been detected in commissioning data of the Sloan Digital Sky Survey (SDSS), providing deep 
five-color photometry in a 2.5° wide band along the equator. The stars in the tails make up a 
substantial part (~ 1/3) of the current total population of cluster stars in the magnitude interval 
19.5 < i* < 22.0. This reveals that the cluster is subject to heavy mass loss. The orientation of 
the tails provides an important key for the determination of the cluster's Galactic orbit. 



Subject headings: globular clusters: individual (Pal5) 
kinematics and dynamics 

1. Introduction 

Globular clusters are self-gravitating stellar sys- 
tems that experience a time- varying tidal potential 
as they orbit through their parent galaxy. Their 
dynamical evolution is driven by internal effects 
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such as stellar evolution, two-body relaxation and 
binary heating, and by external effects induced by 
the galactic force field, i.e., heating by tidal shocks 
during disk and bulge passages and tidal stripping. 
Both internal and external effects should lead to 
a permanent loss of cluster members and to the 
eventual dissolution of the cluster. The Galactic 
globular clusters observed today are believed to be 
survivors from an initially much more numerous 
population. They appear to be in various stages 
of evolution and dissolution, depending on their 
initial conditions and their galactic orbits (e.g., 
Chernoff & Weinberg 1990, Djorgovsky & Meylan 
1994). Numerical simulations predict that possi- 
bly as many as half of the present-day Galactic 
globulars will not survive for another Hubble time 
(Gnedin & Ostriker 1997). Observational confir- 
mation of the gradual dissolution of globular clus- 
ters and determination of their mass loss rates is 
important in itself, but can also shed light on the 
formation history and structure of the Galactic 
halo and provide constraints on the Galactic po- 
tential. 
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Table 1 

Fundamental parameters of Pal 5. 



Parameter 


Value 


Reference 


a, 8 (J2000) 


229.02°, -0.11° 




l,b 


0.9°, +45.9° 




d 


23.2 kpc 


Harris (1996) 


M 


1.3 ■ 1O 4 M 


Sandage & Hartwick (1977) 


r c 


2.9', 20 pc 


Trager et al. (1995) 


v r 


—56 kms~ 1 


Smith (1985) 


fia cos <5, fi$ 


(-2.44, -0.87) masyr -1 


Schweitzer ct al. (1993) 


fl(X COS <5, /X5 


(-1.0, -2.7) masyr" 1 


Scholz et al. (1998) 


He cos <5, n s 


(-2.55,-1.93) masyr- 1 


Cudworth (priv. comm.) 



First signs of the existence of tidal debris 
around Galactic globular clusters were found in 
studies of globular cluster radial surface density- 
profiles. In many cases the observed profiles de- 
viate from the profile of a best-fit King model at 
the outermost radii and extend beyond the con- 
ventional limiting radius set by the model (e.g., 
Grillmair et al. 1995; Lchmann & Scholz 1997). 
Since the debris produced by a dissolving globu- 
lar cluster will have similar orbital parameters as 
the cluster itself it is expected to drift away from 
the cluster in a leading and a trailing tail which 
are aligned with the cluster's orbit. Searches for 
such tidal tails in the vicinity of Galactic globu- 
lar clusters require large area coverage and have 
so far been conducted with photographic mate- 
rial (Grillmair ct al. 1995; Leon et al. 2000; Testa 
ct al. 2000). While spatially distinct star count 
overdensitics were detected around many clus- 
ters, density fluctuations caused by distant galaxy 
clusters, variable foreground reddening or pho- 
tographic plate inhomogeneities have been found 
to be serious contaminants and have left the lo- 
cation and shape of the tidal tails in most cases 
uncertain. 

The Sloan Digital Sky Survey (SDSS; sec York 
et al. (2000); Gunn ct al. (1998); Fukugita et 
al. (1996)) is designed to provide homogeneous, 
deep (R ~ 23) CCD imaging in five passbands 
(u 1 , g', r',i', z') with contiguous coverage of 10 4 D° 
in the Northern Galactic Cap. Such a high-quality 
database is ideally suited for studies of Galactic 
structure (e.g., Yanny et al. 2000; Ivczic ct al. 
2000) and lends itself in particular to searches for 
tidal debris with unprecedented depth and spatial 



coverage. 

The SDSS commissioning scans happen to cover 
Palomar 5, a remote, sparse halo cluster. The pe- 
culiar characteristics of this cluster, i.e., very low 
mass, large core radius, and low central concentra- 
tion (see Tab. 1) suggest that Pal 5 may have lost a 
significant amount of its initial mass and perhaps 
be close to disruption. Its radial velocity combined 
with determinations of its absolute proper motion 
suggest an eccentric orbit with passages through 
the Galactic disk inside the solar circle (Scholz et 
al. 1998). 

In view of these indications we have searched 
the SDSS data for tidal debris around Pal 5. In 
Section 2 we describe our search method. Results 
are presented in Section 3. 

2. Search method 

In the current phase of SDSS, the region of Pal 5 
is covered by an equatorial band with declinations 
— 1.25 < 8 < +1.25. Our investigation concen- 
trates on the range 226° < a < 232°, i.e., 3° 
east and west of the cluster, and on objects clas- 
sified as SDSS point sources. The typical accu- 
racy of the relative photometry in this data set 
is better than 0.03 mag for stars brighter than 
20.0 mag in g* , r* , or i* , and is ~ 0.1 mag at 
u* = 21.0,3* = 22.2, r* = 21.8, i* = 21.4, and 
z* = 20.0 mag (* is used to indicate that the pho- 
tometric zcropoints are still preliminary). In or- 
der to enable corrections for interstellar extinction 
the SDSS data provide individual extinction val- 
ues based on the maps of Schlegel et al. (1998). 
These values (which are small in the field of Pal5) 
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Fig. 1. — Color-magnitude diagrams for stellar SDSS sources in the field of Pal 5. From left to right: 1. 
(c2, i*) for objects within 3' of the cluster center (core). 2. (ci, i*) for objects within 3' of the cluster center 
(core). 3. (ci,i*) for objects in two circles centered on the cluster tails 4. (c\,i*) for objects in two circles 
in the outer field Note that the field sizes of the tail and field samples shown here are 8 x that of the core 
sample. The thick lines indicate our photometric selection criteria for cluster member candidates, the thin 
lines give 2<7-limits for the expected dispersion of cluster stars in c\ and c 2 due to photometric errors. 



have been subtracted from the observed stellar 
magnitudes. 

One straightforward method for tracing the 
spatial distribution of cluster stars in a large field 
around a cluster is as follows: (i) Define the 
distribution of the cluster population in color- 
magnitude space by means of the stars near the 
cluster center, and that of the field stars by stars 
from a sufficiently distant part of the field, (ii) 
Determine the significance of each point in color- 
magnitude space by means of the local signal-to- 
noise ratio s of the expected true number of cluster 
stars, (i.e. the counted number of cluster stars mi- 
nus the estimated number of contaminating field 
stars, see eq.(2)) in the surroundings of this point, 
(iii) Select those regions in color-magnitude space 
with high significance, i.e., with local signal-to- 
noise s above a certain threshold, extract from the 
entire sample all stars which lie in these particular 
color-magnitude regions, and analyze the spatial 
configuration of these stars. 

This way of empirical photometric filtering was 
first described and used in the two-color study 
by Grillmair et al. (1995). ft could, in princi- 
ple, be easily generalized to a higher dimensional 
photometric space and thus be directly applied 
to the multimagnitude or multicolor-magnitude 
data of SDSS. However, in the case of the very 



sparse cluster Pal 5 it turned out that the num- 
ber of cluster stars was too small to provide rea- 
sonable star count statistics in multidimensional 
color-magnitude cells. In order to avoid problems 
with noise from small numbers we thus worked 
with projections on color-magnitude planes, viz. 
those spanned by the magnitude i* and the color 
indices u* - g*,g* - r*,r* - i*, and i* - z* or 
orthogonal combinations of the latter. 

From g* — r* and r* — i* we tailored new color 
indices ci, C2 as given in eq. (lab), with the advan- 
tage that for cluster stars the systematic variation 
of color with magnitude is mostly contained in C\ 
while variations in C2 are mostly random, i.e. due 
to natural spread and observational errors (see left 
panels of Fig.l). 

ci = 0.907(g*-r*) + 0.421(r*-i*) (la) 
c 2 = -0.421(3* -r*) + 0.907(r* - f) (lb) 

The indices u* — g* and i* — z* were considered 
separately because they are less accurate than the 
others due to decreased sensitivity in the bluest 
and reddest photometric band. 

As the first step of photometric filtering we ap- 
plied simple magnitude-dependent color cuts in 
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Fig. 2. — Projected spatial distribution of the photometric cluster member candidates selected with the 
method described in Section 2. Left: individual star positions. Right: contours of the surface density 
derived by adaptive kernel estimation. Contour levels are 0.19, 0.21, 0.24, 0.3, 0.4, 0.6, 1.0, 2.0, 4.0 and 8.0 
stars/D'. The background has a mean of 0.13 stars/d'. 



u* — g* , i* — z* , and ci according to the borders 
of the cluster population defined by stars from the 
cluster core (r < 3'), and according to the esti- 
mated spread due to photometric errors (for ci cf. 
leftmost panel of Fig. 1). 

After this preselection, we used Grillmair's 
method for filtering in the color-magnitude plane 
of (ci,i*). The local signal-to-noise s of the ex- 
pected true number of clusters stars was deter- 
mined on a grid of mesh size 0.01 mag in ci and 
0.05 mag in i*. For each grid point (index k) 
the number n c (k) of stars in a circle of radius 12' 
around the cluster center and the number nj(k) 
of stars at more than 2° angular distance from the 
cluster center were counted in a color-magnitude 
box of width 0.09 mag in c\ and 0.35 mag in i* 
centered on that point. With n c representing the 
number of cluster stars plus underlying field stars, 
rif the number of field stars, and q the ratio of the 
areas on which nj and n c have been sampled, s 
was calculated as given by eq.(2). 

.(*) . y [k i::y k L m 

^n c {k) + q z n f (k) 

The size of the color-magnitude boxes and the 
overlap between boxes around neighboring grid 
points assure that s is a sufficiently smooth func- 
tion. From the local signal-to-noise s one obtains 
a filtering mask in the (ci,i*) plane by setting a 



threshold su m < s max and by isolating the re- 
gion in the grid (around the maximum of s) with 
s > sum- In order to find the most appropriate 
mask, we went through a scries of gradually de- 
creasing thresholds, counted for each threshold the 
cumulative number of stars in the corresponding 
mask in the area of the cluster's tails (N t ) and in 
the outer field (Nf), and determined from these 
numbers the cumulative signal-to-noise ratio SNR 
of the expected true number of cluster stars in the 
area of the tails (eq.(3)). 

SNR = ^3 (3) 

y/Nt+W~ 2 Nf 

The filtering mask was then chosen such that 
the cumulative signal-to noise reaches a maximum. 
As shown in Fig. 1 (second panel from left) this 
mask cuts out the zone from the bottom of the 
subgiant branch to the main-sequence turnoff and 
further down the main sequence to i* = 22.0 mag. 
In the range 19.5 < i* < 21.5 the width of the 
mask approximately coincides with the 2a limits 
for the dispersion of cluster stars in c\ as derived 
from the median values of the estimated photo- 
metric errors. The two panels on the right in Fig. 1 
give an example of the detection of cluster mem- 
ber candidates outside the cluster using the filter 
mask in the area of the cluster's tails and in the 
area of the outer field. 
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The spatial configuration of the complete sam- 
ple of member candidates obtained by the pho- 
tometric filtering process is shown in Fig. 2 (left 
panel). In the final step of data processing, the 
distribution of individual star positions was trans- 
formed into a smooth surface density function by 
means of adaptive kernel estimation (Silverman 
1986). A standard parabolic kernel was used, with 
the kernel radius set to the angular distance of the 
70th nearest neighbor of each star. This yields the 
surface density distribution shown in the contour 
plot of Fig. 2 (right panel). 

3. Discussion 

3.1. Characteristics of the tails 

Fig. 2 shows that the density enhancements of 
point sources around the cluster form two spectac- 
ular tails which emerge from the cluster in north- 
ern and southern direction and turn over to the 
north-east and south-west, respectively, at angu- 
lar distances of ~ 0.2° (80 pc in projected linear 
distance) from the cluster center. The tails stretch 
out almost symmetrically to both sides and exhibit 
clumps at a distance of ~ 0.8°, i.e., 320 pc from 
the cluster. In total, the tails are visible along an 
arc of 2.6°. A weaker clump at the southern edge 
of the current field suggest that the tails might in 
fact continue to even larger distances. 

In the two big clumps, the surface density of 
stars that fall inside our colour magnitude filter 
is about 2.3 times as high as in the surround- 
ing field. Summing up the number of stars above 
background in the region of the tails and com- 
paring them to the stars within a circle of radius 
r < 12' around the cluster center we find that 
within our color-magnitude window the tails con- 
tain ~ 0.48 times the number of stars in the clus- 
ter. In other words, the tails comprise ~ 32% 
of the currently detected total population of clus- 
ter stars at and below the main-sequence turnoff. 
This is a rough (but conservative) estimate be- 
cause the object is seen in a non-face-on projection 
which does not reveal a clear border between clus- 
ter and tail. Nonetheless it gives impressive evi- 
dence for heavy mass loss, confirming conjectures 
drawn from the low mass and low concentration 
of the cluster. 

The structure of the observed tails follows the 
principal expectations for tidal tails and closely 



agrees with the results of recent N-body simula- 
tions of globular cluster tides (Combes et al. 2000) . 
Basically, cluster members drift to the outer part 
of the cluster after acceleration in disk or bulge 
shocks and leave the cluster in the vicinity of the 
(Lagrange-) points of force balance between the 
cluster and the tidal field, i.e., in the direction to 
the galactic center and anticenter. Due to differen- 
tial galactic rotation their trajectories then bend 
around and continue approximately parallel to the 
orbit of the cluster. The appearance of clumps in 
the tails is also supported by numerical simula- 
tions. They can either be caused by the enhanced 
release of particles after strong shock events or by 
caustics of the trajectories in phase space. 

3.2. The role of contaminants 

In view of the striking resemblance of the de- 
tected structures to the expected properties of 
tidal tails, significant contamination by clustered 
background objects seems a priori unlikely. Nev- 
ertheless we checked this point by analysing the 
density and colors of non-pointlike SDSS sources 
around Pal 5 which by their shape are classified 
as galaxies. Their spatial distribution is clumped 
and reveals known galaxy clusters like Abell 2050 
and 2035. However, most of these sources do not 
fall into the color-magnitude window for members 
of Pal 5. If our selection criteria for Pal 5 mem- 
bers are applied to the galaxy sample, the surface 
density of galaxies drops to the level of 0.3 to 0.6 
of the field background density of the stellar sam- 
ple. Moreover, the pattern of density variations 
in the galaxy sample does not correlate with the 
location of the tidal tails. Therefore, objects like 
those in the galaxy sample are not likely to cause 
significant disturbances in the stellar sample. The 
only remaining contaminants are compact galaxies 
with bluer colors that may not be well represented 
in the sample of known galaxies. We believe that 
such objects have mostly been eliminated by our 
color cut u* ~ g* < 0.4 mag. Finally, fluctuations 
in stellar surface density due to variable interstel- 
lar absorption can be ruled out because the mean 
level of absorption in the region around Pal 5 is 
low and there is no hint for strong variations (val- 
ues of Eb-v m the maps of Schlegel et al. (1998) 
range between 0.05 mag and 0.07 mag). 
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Fig. 3. — Contour plot of the surface den- 
sity of cluster candidates in galactic coordinates 
(I cos b, b) overlayed with different orbital paths of 
the cluster according to different determinations of 
its absolute proper motion: S93 = Schweitzer et 
al. (1993); S98 = Scholz et al. (1998); C98 = Cud- 
worth (unpublished), cited in Scholz et al. (1998). 
The solid line presents our improved estimate of 
the orbit based on the geometry and orientation 
of the tidal tails (fixing the direction of tangential 
motion) and the proper motions by Cudworth and 
by Scholz et al. (estimating the tangential veloc- 
ity)- 



3.3. Implications for and from the clus- 
ter's orbit 

The orientation of the tails provides unique in- 
formation on the direction of the cluster's tan- 
gential motion since it is known that the lead- 
ing and trailing parts must fit in with the inner 
and outer side of the local orbit, respectively. Fig. 
3 reveals that the tails stretch out in the direc- 
tion of constant b, and that the tangential mo- 
tion is very likely westward (prograde rotation). 
The absolute proper motions by Schweitzer et al. 
(1993) and Scholz et al. (1998) yield very differ- 
ent predictions for the local orbit, although at 
least the sense of rotation is in agreement. The 
proper motion given by Cudworth (priv. comm., 
unpublished revision of the work by Schweitzer 
et al.) yields a local orbit which lies closer to 
the tails. From the observed orientation of the 
tails we estimate that the tangential velocity vec- 



tor points ~ 15° north of the line of constant 
b. In order to meet this constraint the proper 
motion needs to be modified by not more than 
0.4 masyr -1 w.r.t. Cudworth's values. We thus 
adopt [ncosb,fib = —0.93, +0.25 mas yr _1 for the 
cluster's proper motion in the galactic rest frame. 
With these values we obtain an orbit (using the 
galactic potential of Allen & Santillan 1991) with 
apo- and perigalactic distances of 19.0 kpc and 
7.0 kpc, and with disk passages at —137 Myr, 
-292 Myr and -472 Myr, taking place at Galacto- 
centric radii of 9.4, 18.4 and, 8.3 kpc, respectively. 
Similar orbits are obtained with the more detailed 
galactic potentials of Dchnen & Binney (1998). 
We tend to believe that the observed overdensities 
close to the cluster result from the latest disk pas- 
sage while the clumps in the tails at distances of 
<~ 0.8° from the cluster might be associated with 
the earlier passage through the inner disk about 
470 Myr ago. This, however, has to be investi- 
gated more thoroughly with N-body simulations. 
The next passage through the galactic disk pre- 
dicted by our model orbit will be in 113 Myr and 
will happen close to perigalacticon. If true, this 
will again produce a strong tidal shock that may 
eventually dissolve the cluster completely. 

3.4. Outlook 

The SDSS will eventually cover a much larger 
region around Pal 5 than currently available. 
Larger area coverage will enable us to constrain 
the orbit and mass loss more tightly. We can fur- 
ther constrain our model orbit by obtaining radial 
velocities of stars in the tails. We predict a local 
radial velocity gradient of 5.7 kms _1 deg _1 , i.e. 
<~ 9 kms -1 difference between the radial velocities 
of stars in the two tidal tail clumps. Since Pal 5 
contains very few luminous red giants, even fewer 
are expected in its tails, thus kinematic studies 
will have to concentrate on fainter stars requiring 
large telescopes. 

The Sloan Digital Sky Survey (SDSS) is a joint 
project of The University of Chicago, Fcrmilab, 
the Institute for Advanced Study, the Japan Par- 
ticipation Group, The Johns Hopkins University, 
the Max Planck Institute for Astronomy, New 
Mexico State University, Princeton University, the 
United States Naval Observatory, and the Univer- 
sity of Washington. Apache Point Observatory, 
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site of the SDSS, is operated by the Astrophysi- 
cal Research Consortium. Funding for the project 
has been provided by the Alfred P. Sloan Founda- 
tion, the SDSS member institutions, the National 
Aeronautics and Space Administration, the Na- 
tional Science Foundation, the U.S. Department 
of Energy, and Monbusho. The SDSS Web site is 
http://www.sdss.org/. 
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